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Abstract. We consider the linearized 2D inviscid shallow water equations in a rectangle. 
A set of boundary conditions is proposed which make these equations well-posed. Several 
different cases occur depending on the relative values of the reference velocities (uo,Vo) 
and reference height 0o ( sub- or super-critical flow at each part of the boundary) . 



1. Introduction 

Motivated by the study of the inviscid linearized primitive equation (as a step towards the 
study of the full nonlinear equations) , we consider in this article the linearized fully inviscid 
shallow water equations in a rectangle. To the best of our knowledge, the well-posedness of 
the linearized inviscid shallow water equations in a rectangle has not yet been addressed. In 
particular, the existence of corners in the geometrical domain is a subject of mathematical 
concern since the seminal works |Osh73l IOsh74| which show the possible occurrence of major 
singularities in the corners for certain choices of the boundary conditions. The boundary 
conditions that we choose do not lead to singularities and allow us to develop results of 
existence and uniqueness of a solution. As explained below, in the rationale of the Local 
Area Models (LAMs) to which this work relates (see e.g. |RTT08al IRTT08bl IWPT97j ). 
we are not looking for all the possible boundary conditions which make the equations well- 
posed, but only certain sets of such boundary conditions, any set in fact. 

In this article, the well-posedness of the linearized inviscid shallow water equations is 
classically conducted using the linear semigroup theory. We write these equations as an 
abstract evolution equation in a Hilbert space: 

with initial data U (0) = Uq given, which means that we prescribe the boundary conditions 
at t = (see below). In this approach, the main difficulty is to show that (AU,U) > 
for all U in the domain V(A) of A (and similarly (A*U, U) > 0, VU G T>(A*)), where A is 
the abstract associated operator acting in L 2 (f2) 3 . This result is achieved using a suitable 
approximation technique. Five different cases occur depending on whether the linearized 
base flow with velocities (uo,vq) and height cpo is sub- or super-critical on each side of 
the rectangle (wq,Vq > or < g<po). The fifth case corresponding to Uq + Vq < g<po is fully 
subcritical (uq,Vq < gcfio); furthermore the time independent part of the shallow water 
equations is partly elliptic in this case, but of course the full shallow water system remains 
hyperbolic. Boundary value problems of the Dirichlet or Neumann type appear in this case. 
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Of course the presence of the corners also raises difficulties for these elliptic problems, but 
the desired regularity of the solutions of the corresponding Dirichlet-Neumann problem is 
achieved using results from |Gri85j . 

The study of the Local Area Models (LAMs) in the atmosphere and oceans sciences, leads 
to initial and boundary value problems for the inviscid primitive equations in these domains, 
and we know since [OS78j (see also |TT03| ) that these problems are complicated leading to 
nonlocal boundary value problems. As explained in e.g. |WPT97] the choice of the boundary 
condition is important for the numerical simulations, as one wishes boundary conditions 
leading to well-posed problems (to avoid numerical explosion) and boundary conditions 
that are transparent, letting the waves move freely inside and outside the domain. In the 
perspective of a LAM imbedded in a large domain, any such set of conditions is acceptable. 

The linearized shallow water equations that we consider read 



(1.1) 



( u t + u u x + v u y + gcj) x - fv = 0, 
vt + u v x + v v y + g<j) y + fu = 0, 
<h + u o<t>x + v (f>y + <f>o(u x + v y ) = 0; 



here U = (u, v, (/))*, (x, y) 6 := (0, L{) x (0, L2), uo,vo,(f>o are positive constants, and g 
is the gravitational acceleration, / is the Coriolis parameter. As explained in e.g. [OS78, 
ITT03| IRTT08b| . the inviscid primitive equations in a cube can be expanded in vertical 
modes expansion producing a set of coupled bi-dimensional equations, similar to coupled 
shallow water equations. For one linearized mode, the system that we obtain is very similar 



to ( 1.1 ) as we show at the end of this article. With this respect the results that we obtain for 
the one-mode linearized primitive equations generalize in part the results from [RTT08b~] 
in which we assumed vq = 0; the full coupled system as considered in [RTT08b] will be 
studied elsewhere. 

As we know the literature on the shallow water equations is very vast, both on the 
theoretical and computational aspects, considering the viscous equations or the partly or 
totally inviscid equations and considering that the height is either always strictly positive or 
that it can vanish. See e.g. [ABBKPl IAB031 IABPSI ISLTTl IBPSTT] on the computational 
side and see e.g. [Ure95l IPT101 IPTlll IHPTlll IRTT08a] on the theoretical side. But as we 
said we did not find this problem studied in the literature. The general hyperbolic equations 
in a domain with corners are also considered in [Tan78, KT80J, but the boundary conditions 
they imposed is different from ours, and the space they used is not the usual L 2 space which 
we use in this article. 

This article is organized as follows. After this introductory section, we describe in Section 
[2] the boundary conditions which one can associate to the time independent part A of the 
shallow water equations (dropping also the Coriolis terms). We also present a first useful 
trace theorem. Then in Section [3] we derive various density theorems, density of certain 
smooth functions in certain function spaces. Section HI is devoted to the fully supercritical 
case when Uq and Vq are both larger than gcf)Q ; Sectionpl is devoted to the mixed cases when 
Uq > gcj)o and v 2 < g<fto or vice versa; Section [6] tackles the fully hyperbolic subcritical case 
corresponding to Uq < gcfto and Vq < gcfto whereas Uq + v 2 > g<j)Q. Section [7] is a technical 
section giving a regularity result for an associated elliptic problem. It prepares to Section 
[8] in which we deal with the mixed subcritical case, where appears a partly elliptic and 
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partly hyperbolic operator A; this case corresponds to Uq + Vq < g4>o, so that Uq < g4>o and 
Vq < g<j)Q. In Section [9] we state the general existence and uniqueness theorem proved by 
application of the Hille-Phillips-Yoshida theorem, which covers all cases. We finally make 
the connection with the equations for one mode of the linearized primitive equations and 
interpret our earlier result for that case. 

We conclude this introduction with a last remark which will be slightly re- 
worded in the final version of this article. This work is at the interface between 
hyperbolic equations and parabolic-elliptic equations; the equations are hyper- 
bolic, the methods used are rather those of parabolic and elliptic equations. We 
believe that the results presented here are new; if they are known, they were 
at least not known by us as we said nor by those working in geophysical fluid 
mechanics who need them most |WPT97j . The authors are not specialists of 
hyperbolic equations and they may have missed some relevant references. They 
will gratefully receive any suggestion from the Editors or from the Referees on 
the bibliography or else. In some sense, this article is a continuation of [RTT08b] 
published in the Journal de Mathematiques Pures et Appliquees which dealt 
with the case where vq = 0. However, the methods here are different. Further- 
more beside the four cases that one can surmise by comparison with [RTT08b], 
there is a fifth unexpected case which is the object of a specific study in Sec- 
tion [8| Finally in |RTT08b|, we related the system to the primitive equations of 
the atmosphere and the oceans. Here we rather relate it to the shallow water 
equations, but the systems are essentially the same, the relation between these 
two systems is clarified in Remark |9.2| , 

2. Boundary conditions 



In order to derive the boundary conditions for (1.1), we introduce the time independent 
2D shallow water equations operator: 

u u x + v u y + g<f> x 

(2.1) AU:=A\v\ = \ u v x + v v y + gcpy 

Uo4> x + V (j) y + (j) (u x + v y ) 





we set 

£i=IO n j , £ 2 

\(f>o n / 

and we aim to determine boundary conditions which are suitable for the system 
(2.2) AU = £iU x + £ 2 U y = F. 

In the following, we assume that no > 0, vo > 0. The case where no and, or no are negative 
can be treated in a similar manner; we do not study the non-generic case where one or 
both of no, no vanish. By its physical meaning (height), 0o > and we do not study the 
case where </>o = 0. Similary, we only study the generic cases where Uq ^ gcfto, Vq 7^ gcj)Q 
and Uq + Vq 7^ g4>Q . As we will see below, we have five cases to study, for four of them 



equations (2.2) are fully hyperbolic and for the fifth one the equations are partly hyperbolic 



and partly elliptic. 
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We first consider the situation when Uq + v g — gcfto is positive, which produces four fully 
hyperbolic cases, and we call them the supercritical case, the mixed case (two sub-cases), 
and the fully hyperbolic subcritical case. Let us set kq = ^ g(u^ + v q — g(po)/<Po- 

2.1. The supercritical case. The supercritical (fully supersonic) case corresponds to the 
case where: 



(2.3) u > g<p , v Q > gcpQ. 

We compute the eigenvalues of £\ and £2- 

Eigenvalues of £\ : uq — \/g4>o, uq , uq + 
Eigenvalues of £2 : — Vd&O' v > v + Vw) • 



(2.4) 



With assumption (2.3), we see that all the eigenvalues of £\ and £2 are positive. Thus, it 
is natural to treat either the x- or y-direction as the time-like direction. Let us choose the 
^/-direction, which means that we first need to specify the boundary conditions at y = 
(time-like initial conditions). Multiplying both sides of (2.2) by £ 2 ~ 1 gives 

(2.5) U y + £ 2 ~ l £iU x = £ 2 X F. 

Directly or on using Matlab, we compute 



£ 2 1 £\ 



v 



g<t>o) 



vo4>o 





u v 

-U (j)Q 



(2.6) 
and write 

(2.7) P- 1 -£ 2 X £ V P = diagfa, A 2 , A 3 ) 

where P has a complicated expression, whereas 

/ 



9 -{vl-9<t>of 



'■(I 



-guo 

UqVq 



P- 



V 




and 
(2.8) 



^0^0 + ^0^0 



A, 



Uo 

vq' 



Al — — I f— A 2 — 2" 

Vq - 9<Po Vq 

Then the eigenvalues Ai,A2,A3 of £ 2 ~ 1 £i are all positive under assumption (2.3). Now we 
introduce the new variables 



' 2k. 



(2.9) 



P- 



/ (V U - U Q V + K (, 
(VQU - U V - K (, 

\[u u + v Q v + g4>)^rj^2 



2fto 



multiplying both sides of (2.5) by P 1 
(2.10) 




P- X £ 2 X F, 
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which is a fully decoupled system. Therefore, from the general Kreiss-Lopatinskii theory 
for hyperbolic boundary value problems (see e.g. |Kre7CH |Lo~p 70 or Chapter 4 in |BS07j ). 



in order to solve the fully decoupled system (2.10) in (£,77, £), it is necessary and sufficient 
to specify the boundary conditions for (£, 77, £) at x = since Ai, A2 and A3 are all positive. 
That is equivalent to specifying the boundary conditions for U = (u, v, (ft) at x = 0. 

In conclusion, we will choose to specify the boundary conditions for U = (u, v, (ft) at x = 
and y = in the supercritical case corresponding to assumption (12 .31) . 



2.2. The mixed case. There are two sub-cases in the mixed case, corresponding to 



(2.11) 

and 

(2.12) 



ul < 



Mn > 



^0 > #o- 



^0, v < g(ft . 



These two cases are similar, we will only consider the first one. 



Under assumption (2.11), the calculations (2.4)-(2.10) are still valid, but £2 only has all 



its eigenvalues positive. Hence, we still choose the y-direction as the time-like direction, and 
thus we need to specify the boundary conditions at y = (time-like initial conditions), as 



in the supercritical case. Of course, we also have (2.10), but the difference here is that only 
Ai,A3 are positive while A2 is negative under assumption ( |2.11 ). Hence in order to solve 
system (2.10) in (£, 77, £), it is necessary and sufficient to specify the boundary conditions 
for £ and £ at x = since Ai and A3 are positive and the boundary condition for 77 at x = L\ 
since Ao is neg ative (see e.g. |BS07| ). 

In conclusion, we choose to specify the boundary conditions for U = (u, v, (ft) at y = 0, 
and the boundary conditions for £ ~ vqu — uqv + KQ(ft and £ ~ uqu + vqv + g(ft at x = 0, and 
a boundary condition for 77 ~ vqu — uqv — kq4> at x = L\ under assumption ( |2.11| ). 



For the other case, corresponding to (2.12), we choose, with similar arguments, to specify 
the boundary conditions for U = (u,v,(ft) at x = 0, the boundary conditions for 77 ~ vqu — 
uov — KQ(ft and £ ~ UQU+VQV+g(ft at y = 0, and a boundary condition for £ ~ v qu — uqv + kq4> 
at y = L2. 



2.3. The fully hyperbolic subcritical case. 

sponds to the case where 



The fully hyperbolic subcritical case corre- 



(2.13) 



^0, v < g(ft . 



> 0, the system (2.5) is still fully 



Since we are still under the assumption Uq + Vq — , 
hyperbolic. But now Ai,A2 < and A3 > 0. Therefore, by analogy with the mixed case 
and due to the energy estimates justified in Section [6j we choose to specify the boundary 
conditions for £ ~ vqu — uqv + KQ(ft and £ ~ uqu + vqv + g(ft at x = 0, for 77 ~ vqu — uqv — Kq4> 
at x = Li, and the boundary conditions for 77 ~ vqu — uqv — KQ(ft and £ ~ uqu + vqv + g(ft 
at y = 0, and for £ ~ vqu — uqV + KQ(ft at y = L2. 



We now consider the situation when Uq 



g(fto is negative, which we call the mixed 



subcritical case. We set K\ = \J g(g<fto — Uq — Vq)/^. 
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2.4. The mixed subcritical case. 

where 

(2.14) 

which implies that 
(2.15) 



The mixed subcritical case corresponds to the case 



Uq + Vq < gfo, 



Uq < 900, vl < gcpo- 



We notice that (2.7) still holds while P is a complex matrix and Ai,A2 are complex num- 



bers, which indicates that there are some elliptic modes hidden in (2.2). In order to avoid 
computations with complex numbers and to take advantage of the regularity results valid 
for elliptic systems (see Sec. [7|), we now define P,t;,T],( by setting 





P 



-i 







v g 

we also denote So = diag(l, 1, g/cpo). 
(ul + ug)P*5 gives 







l — l 






) 1 


V 


V 


H 













VqU — UqV 
Kl<j) 
UqU + V V + ( 



Then multiplying to the left both sides of (2.2) by 





( Uq 


gyp 





(2.16) 


gyp 


-Uq 







V o 





Uq 





(ul + vDP'SoF. 



We see from (2.16) that the equation satisfied by £ is decoupled, and is a transport 
(hyperbolic) equation. Therefore in order to solve for £, we need to assign the boundary 
conditions for ( at x = and y = since uq and vq are both positive. 



From (|2.16|) again, the equations satisfied by £ and r\ read 



(2.17) 



Uq 

gyp 
«i 



gyp 

-Uq 



+ 



gup 

Kl 



gyp 

Kl 



where /i,/2 are the first two components of (uq + vl)P t SoF. The equations (2.17) are 
elliptic since the four constants uq,vq, gvo/nx, —guo/ni satisfy condition ( |7.2[ ) (see Sec. [7]). 
Hence, by Theorem 7.1 and Remark 7.1, in order to solve (2.17) for £ and 77, we can assign 



different kinds of boundary conditions for £ and rj as long as they satisfy condition (7.1). 



However, as we will see in Section [8j by taking into account the energy estimates, only 
some of the boundary conditions satisfying (7.1 ) are suitable for the well-posedness. In this 
case we will assign the boundary conditions for £ at x = and y = 0, and the boundary 
conditions for r\ at x = L\ and y = L2. 

In conclusion, we choose to specify the boundary conditions for £ = vqu — uqv and 
£ = uoii + vqv + 50 at x = and y = 0, and the boundary conditions for ~ u at a; = L\ 
and y 



1/2 when assumption (2.14) holds. 



2.5. A trace theorem. Since we have to assign boundary conditions on 80,, we first need 
to make sure that the desired traces at the boundary make sense. Thus we consider the 
space: 

X(Q) = {U £H = L 2 (Q) 3 ,AU € H}, 



endowed with its natural Hilbert norm 
trace result. 



\U\\h + 



\AU\\ 2 H )Vi. 



Then we have the following 
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Proposition 2.1. If U = (u,v,<f)) t G X(Q), then the traces of U are defined on all of 
dVL, i.e. the traces of U are defined at x = 0,L\, and y = 0,L 2 , and they belong to the 
respective spaces H~ 1 (0,L2) and H^ 1 ^, L\) . Furthermore the trace operators are linear 
continuous in the corresponding spaces, e.g., U G X(£l) — > U\ x= q is continuous from X(il.) 
into H- 1 (0,L 2 ) 3 . 



Proof. Since U G L 2 (Q) 3 = L^(0, L\\ L 2 (0, L 2 ) 3 ), we see that U y 
L 2 X (0, Li; Hy 1 ^, L 2 f). Since AU G L 2 (n) 3 , that is 



dU/dy belongs to 



(2.18) 



u u x + v Uy + g4> x 
uqVx + VQVy + g(f) y 

Uq4>X + VQCjiy + (f)0(u x + Vy) 



GL2(0,L i; L2(0,L 2 )), 
eL 2 x (0,L 1 ;Ll(0,L 2 )), 
eLl(0,L i; L 2 y (0,L 2 )), 



we infer from (2. 18) ^3 that 
(2.19) 



u u x +g(f) x G L 2 (0, L\\ H 1 (0, L 2 )), 



UqPx + 



Noticing that the determinant of the coefficient matrix in (2.19) is nonzero since we exclude 
the non-generic cases where Uq = gcfro, we obtain that u x ,cf) x G L^.(0,Li;H~ (0,L2)). We 
also find that v x G l£(0, L x ; -ff^ x (0, L 2 )) from (2.18) 2 . Therefore, we have 

u x ,v x ,4> x G L 2 x {0,Li;H y 1 {Q,L 2 )) 1 

and in combination with u,v,(f) G -Z/^(0, Li; L 2 (0, L 2 )), we obtain that 



(2.20) 



« J ^ J ^GC :c ([0,L 1 ];ff- 1 (0,X 2 )), 
are well-defined at x 



which shows that the traces of u,v,4> are well-defined at x = and Li, and belong to 
iif~ 1 (0, L 2 ). The continuity of the corresponding mappings is easy. The proof for the traces 
at y = and L 2 is similar. □ 



3. The density theorems 

In this section, we establish general density theorems regarding functions defined on the 
domain Q = (0, L\) x (0,L 2 ). This theorem will be needed for proving later on that —A is 
the infinitesimal generator of a semigroup of contraction, where A is the 2D shallow water 
equations operator A associated with the suitable boundary conditions. 

For A fixed, A G M., A 7^ 0, we set TO = 9 y + A^, and introduce the function space 

ATi(fi) = {6 G L 2 (n), T6 = 6 y + \6 X G L 2 (Q)}. 

We need to show that the smooth functions are dense in X\(Q). Later on we will prove 
more involved density theorems, showing that if u G Ai(fJ) vanishes on certain parts of d£l, 
then u can be approximated in X\{yi) by smooth functions, vanishing on the same parts of 
the boundary. For the moment, we prove the following. 

Proposition 3.1. C°°(fT) n Xi{9) is dense in Xi(Q). 
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Proof. Using a proper covering of £7 by sets Oq, 0\, ■ ■ ■ , On, we consider a partition of unity 
subordinated to this covering, 1 = X^^nV'v Here and again in this section we will use a 
covering of ft consisting of Oq, a relatively compact subset of f2, and of sets Oi of one of the 
following types: Oi is a ball centered at one of the corners of Q, which does not intersect 
the two other sides of 0; or Oi is a ball centered on one of the sides of Q which does not 
intersect any of the three other sides of O. 

If 8 G Xi(Q), then Oipi G X\(Q), so that we only need to approximate Oipi by smooth 
functions. Here the support of tpi is contained in the ball Oi, and we start with considering 
the ball Oq, relatively compact in J7, then considering the balls Oi on the boundary d£l. 

For any function v defined in Q, here and again in the following we denote by v the 
function equal to v in £1 and to outside f2. We first consider the case where ipi = ipo, and 
Oi = Oq, which is relatively compact in S7. Let p be a mollifier such that p > 0, f p = 1, 
and p has compact support. For v = Oipo, the standard mollifier theory (see e.g. Appendix 
C in [Eva98j ) shows that for e — > 0, we have 



(3.1) 



p e * V — > V, 

p e * Tv — > Tv, 



in L 2 (M 2 ). Now observing that Tv = Tv, we find 
(3.2) T(p e *v) = p e *Tv = p e *Tv -+Tv, 

in L 2 (M 2 ) as e — > 0. Since Oq is relatively compact in $7, for e s mall enoug h, p € * v is 
supported in Q. If v e = p € * v, then v e \ n converges to v in Xi(Q) by (3.1 )i and (3.2). 



We then consider the case where ipi = ipi, and Oi = 0\ which is a ball centered at the 
origin (0,0); the other cases are similar or simpler. Let p be the same mollifier as before, 
but now p is compactly supported in {x < 0,y < 0}. Then for v = Otpi, we observe that 

(3.3) Tv = Tv + n, 

where p, is a measure supported on {x = 0} U {y = 0}. Then by mollification with this p 
(see |Hor65] ). we observe that, if v e = p e * v, then 

(3.4) Tv e = p e *Tv + p e * p, 

where p e * p is supported in QP by the choice of p. Then, as e — > 0, v e — > v in L 2 (R 2 ), and 



(3.5) 



v e \ n —> v, in L 2 (Q); 
Tv e L^ Tv, in L 2 (ty, 



\Q 

which shows that v e \ Q converges to v in Xi(Q). □ 



Remark 3.1. Proposition \3.1\ is also valid with T replaced by any first order differential 
operator with constant coefficients. 

We are now going to prove the density theorems involving the boundary dVl, and we will 
successively consider the scalar case and the vector case. Here and throughout this article 
we denote by Fi,r2,r3,r4 the boundaries x = 0, x = L\,y = 0,y = Li respectively, and 
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let r be any union of the sets ri,r2,r3,r4. In the scalar case, we introduce the function 
spaces: 

(3.6) Afr(fi) = {0 G L 2 (Q),T6 = 6 y + X9 X G L 2 (Q), 8\ r = 0}, 

Vr(^) = {8 G C°°(r2), and vanishes in a neighborhood of T}. 
We then state the density theorem: 

Theorem 3.1. Suppose that T = T1UT3 and A > 0, then we have 

Vr(^) n Xr(P>) is dense in X r (£l). 



To prove Theorem |3.1[ we consider the same partition of unity as before; we are led to 
approximate functions like 6ipi by functions of Vr(^) H A'r(fi). The proof for OtpQ is easy as 
before, we will then consider other functions 9ipi depending on the type of the support 0j. 

In order to obtain the local density results regarding certain parts of the boundaries, we 
introduce the new coordinate system (x',y') such that 



which is equivalent to 



1 + A 2 



A 1 

-1 A 



We denote by the image of by this transformation for all i G {1, 2, 3, 4}; we also denote 
by Q' the image of Q, and denote by 0' the transform of 8. Let (ay, y') and (/?,/, y') be the 
end points of the intersection of f2' with the line y' fixed, (a y ',y') G 1^ or T' 3 , (f3 y >,y') G T' 2 
or r' 4 (see Figure 3.1). 



1 

\L<2 




-1/ 






/ a v' 


AN 




y' 


j 















Figure 3.1. The domain f2' 



Direct computation shows that 



d 1 , d d v 

dx' 1 + A 2 dy dx 
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Therefore, if 9 G Ar(f2), we see that 9' and 9 ' , both belong to L 2 (Q'), and 9'(a y /,y') = in 
the new coordinate system (x 1 ', y'). Now for fixed y' , we have 

(3.7) 0'OW) = / fAx^dx. 



Raising (3.7) to power 2 and integrating with respect to y' from —L\ to AL2, then using 



the Cauchy-Schwarz inequality, we finally arrive at 

(3-8) / \9>(f3 y ,,y')\ 2 dy'<C\\9> xl \\ 2 L2{n , y 

where C is a constant independent of 9' . Transforming back into the original system (x,y), 

(3.8) is equivalent to 

(3.9) J 1 9{x, L 2 ) 2 dx + £ 2 9(L 1 ,y) 2 dy < C||T0||* 2( 



Inequalities (3.8) and (3.9) say that the traces 9' L, r , and 9 



V v v> A "ir 2 ,r 4 
respect to the corresponding segments, e.g. 0'| r , G i 2 (l"4 



are L 2 functions with 



Now we are ready to show the local density results near the boundaries in the case (3.6) 



There are eight cases to consider, and we prove them in the following lemmas corresponding 
to four balls O centered at the four corners and four balls O centered on each of the sides. 

Lemma 3.1. Suppose 9 = 9(x,y) G A^r(O) and supp9 C0nH, where O is a ball centered 



on T3, which does not intersect Ti U T2 U T4 (see Figure 3.2). 



L 



Li x 

Figure 3.2. Ball centered on T 3 



Then we have 
(3.10) T9 = T9. 

Assume that p G V(M. 2 ) has its support compact in {\x\ < y}; then p € * 9 vanishes in a 
neighborhood of T, and as e — > 0, 

Up e *9)\ g ^9, mL 2 (n), 
\T((p e *9)\ n )^T9, inL 2 (n). 



(3.11) 



Proof. In the new coordinate system (x' ,y'), 9' has its support compact in O' f~l CI', where 
O' , the image of O, is a ball centered on T' 3 , which does not intersect U r' 2 U T' 4 . We 
extend 9' by for x' + Xy' < (y < 0), and denote this extension by 9'. Noticing that 
9' = on T' 3 = {x' + Xy' = 0}, we then have 

(3.12) 9' X ,=Z, 
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since 9 1 , belongs to L 2 (f2'). Transforming back to the original system (x,y), (3.12) implies 
that 

(3.13) TO = TO. 



Now, for e — > 0, with the same arguments as for (3.1) and (3.2), and using (3.13), we find 

( Pe *9^0, _ 
1 T( Pe * 9) -> TO, 



(3.14) 



in L 2 (R 2 ). 

Since the support of p e *9 is contained in supp(p e ) + supp(0) , we infer that supp{p e *9) is at 
a distance away from the x-axis, and included in the interior of $7. Hence, p e *9 vanishes in a 
neighborhood of T by considering the position of O, so that (p e * 9) L E Vr(O); furthermore 



(3.14) implies that 



Up e *§)\ 9 ^e, 

\T((p e *9)\ n )^T6, 



in L 2 (S7) as e — > 0. The proof is complete. 



□ 



Proceeding exactly as for Lemma 3.1 we obtain again (3.11) when O is a ball centered 



on Ti or at (0,0); in these cases we respectively take p supported in {\y\ < x} and in 
{0 < t;X < y < 2x}. We thus omit the proof here. 

Lemma 3.2. Suppose = 9(x,y) £ Xr(Q) and suppO C O n where O is a ball centered 



on T4, which does not intersect T\ U T2 U T3 (see Figure 3.3). 



ii x 

Figure 3.3. Ball centered on L4 



Then 

(3.15) T9 = T9 + n, 

where p. is a measure supported by = {y = L2}. Assume that p G D(M 2 ) 
compact in {\x\ < y}; then p e * 9 vanishes in a neighborhood ofT, and as e 



has its support 
^0, 



(3.16) 



f(Pe*0)| 
\T{(Pe * 



in L 2 (n), 
T9, inL 2 (n). 



Proof. In the new coordinate system {x' , y' 
O' is a ball centered on V 



9' has its support compact in O' (1 !]', where 

We extend 9' by 



4 , which does not intersect U T' 2 U r' 3 . 
for x' + Ay' > (1 + X 2 )L2 ( y > L2), and denote this extension by 9'. Observing that 
9'\ , e L 2 (T' 4 ) (see (E8|) and fcM), and that ' , belongs to L 2 {9!), we thus have 

4 



(3.17) 
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where p! is a measure supported by (= 9'\ T/ 5-p>). Transforming back to the original 



system (x,y), (3.17) implies that 

(3.18) T9 = T9 + p, 

where p is a measure supported by With the choice of p, p e * p is supported outside of 



U. Hence, for e — > 0, using the same argument as for (3.4) and (3.5), and noticing (3.18) 
we obtain 



(3.19) 



in L 2 (n), 
> T9, in L 2 {Q). 



Finally, p e * 6 vanishes in a neighborhood of T since the support of p t * 6 is away from 

r. □ 



Arguing exactly as for Lemma 3.2, we also obtain (3.16) when O is a ball centered on T2 
or at (L\,L2), by taking respectively p supported in {\y\ < x} and in {0 < \x < y < 2x}. 
Thus, we also omit the proof here. 

We now turn to the cases where O is either centered at (0, L2) or {L\, 0), and both cases 
necessitate a specific treatment. 

Lemma 3.3. Suppose 9 = 9(x,y) £ Xr(£l) and supp9 C O HQ, where O is a ball centered 
at (0, L2), which does not intersect T2 U T3 (see Figure 3.4)- 




ii x 

Figure 3.4. Ball centered at (0,L 2 ) 



Then we have 
(3.20) T9 = T9 + p, 

where p is a measure supported by = {y = £2}- Assume that p E D(IR 2 ) has its support 
compact in {0 < \x < y < 2x}; then p € * 9 vanishes in a neighborhood ofT, and as e — > 0, 

f3211 U Pt *9)\ n ^9, mL\Sl), 

1 ' ' \T((p t *9)\ n )^T9, inL 2 (n). 

Proof. In the new coordinate system (x',y'), 9' has its support compact in O' n where 
O' is a ball centered at (L2,AL2); it does not intersect V 2 U r' 3 . We first extend 9' by 
for Ax' — y' < {x < 0), and denote this extension by 9' . Noticing that 9' = on 
= {Ax' — y' = 0}, we see that 9' x / = 9' x , in the distribution sense in VL' = {x' + Xy' < 
(1 + A 2 )L 2 } (we also see that &~ x , G L 2 {d') because 9' x , £ L 2 (Q,')). We then extend 9' by for 
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x' + Xy' > (1 + A 2 )L 2 (y > L 2 ), and denote this extension by 9'. Noticing that 9'\fr G L 2 (V 4 
since 9' L, G -L 2 ^^), where is the whole line {x' + Ay' = (1 + A 2 )L 2 }, we thus have 



(3.22) 



where // is a measure supported by r 4 . The rest of the proof is the same as for Lemma 
we thus omit it here. 



3.2 



□ 

using 



The case where O is a ball centered at (Li, 0) is treated similarly as in Lemma 3.3 
a modifier p supported in {0 < \x < y < 2x}. We omit the proof. 

We finished considering all the local density results. Thus, at this stage, we are able to 
prove Theorem 3.1 

Proof of Theorem 3.1 Recall the covering sets Oq,Oi,--- ,On described before and the 



partition of unity subordinated to this covering 1 = J2iLo rfi introduced in Proposition 
Given 9 G Xr(Cl), we set 9i = ipiO, so that 0i G Afr(fi). For any e > 0, By Lemmas 3.1 
we find 9\ G C°°(U) such that 



3.1 



3.3 



(3.23) 



Il 2 (q) 



<e, \\T9I 



r ^IL 2 (c) - e ' 



and all 9fs vanish in a neig hborhood of V. Set 9 € = Yli=o e !, then ° e belongs to C°°(Q), 
and vanishes in a neighborhood of T, whic h sho ws that 9 e belongs to Vr(fi). 



Noticing that 9 = ^2*L ^9 



L 2 (n) 



^ =0 9i, (|3.23D implies 




— 2^j=o 



L2(n) <(N + l)e, 



- T9\\ L2{n) < Zlo \\T0! - T9 t \\ L2{n) <(N + l) e , 

which shows that 9 e belongs to A'r(O) and converges to 9 in Xr(Q) as e — )• 0. Hence, the 
proof is complete. □ 



Remark 3.2. Looking back carefully at the proof of Theorem 3.1, we see that Theorem 3.1 
is also valid in the following three cases: 

'r = T 2 U T 4 and A > 0, 
r = Ti U T 4 and A < 0, 
r = T 2 U T 3 and A < 0, 



(3.24) 



provided we choose properly the support of the mollifier. This remark will be useful for 
Theorem 3.2 and when we study the adjoint A* of A later on. 



We now turn to generalizing Theorem 3.1 to the vector case. 

Theorem 3.2. Let n be a positive integer, and suppose that A = diag(Xi, ■ ■ ■ ,A n ) is a 
nonsingular diagonal matrix, and Mi,M 2 ,Q are arbitrary nonsingular matrices satisfying 
Q~ l KQ = M^Mx. Set f = (T 1 , • • • , T n ) where 

r = / riLjr3 orr 2ur 4 , if\ >o, 
|r 1 ur 4 or r 2 ur 3 , if\i <o ; 
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and introduce the function spaces in the vector case: 

x f (n) = {e = (e u ■ ■ ■ A) € L 2 (n) n , MxQ x + M 2 e y e L\n) n , ge| f = o}, 

Vp(fi) = {Q G C°°(f2), and QQ vanishes in a neighborhood ofT}. 
Then we have 

Vf(£l) n Afp(O) zs dense in Xp(Q). 

Proof. Setting G = QQ, we first have G^ + AQ X = QM 2 ~ l {Mi<c> x + M 2 Q y ), and then 
applying Theorem 3.1 or Remark 3.2 to each component of Q and transforming back to G, 
we obtain the result for Theorem 13.21 □ 



(4.2) 



4. The supercritical case 

In this section, we consider the stationary 2D shallow water equations operator for the 
supercritical case; thus we assume that 

(4.1) ul > gfo, Vq > g4> . 

Here, we choose the following boundary conditions according to the discussion in Sub- 
section 12.11 

{u = v = 4> = 0, on Ti = {x = 0}, 
u = v = <j> = 0, on T3 = {y = 0}. 

We then define the unbounded operator A on H = L 2 (Q) 3 , with AU = AU, MU £ V(A) 
and 

V(A) = {U G H = L 2 (nf,AU G H, and U satisfies Q}. 
We also introduce the corresponding density boundary conditions: 
(4.3) U vanishes in a neighborhood of T\ UT3, 

and the function space: 

V(ft) = {U G C°°(n) 3 , and U satisfies B}. 



Setting T = (ri U T3, Ti U T3, T\ U 1^3), we have that P 1 U\f = is equivalent to (4.2 ) and 
P~ 1 U vanishing in a neighborhood of T is the same as (4.3), where P is as in Subsection 



2.1 Then we obtain the following result from Theorem E2fwith Q = P' 1 , M x = £1, M 2 = 
£ 2 , A = diag(\i, A 2 , A 3 ). 

Lemma 4.1. V(O) f]V(A) is dense in V(A). 

4.1. Positivity of A and its adjoint A*. We endow the space H = L 2 (0) 3 with the 
Hilbert scalar product and norm: 

9_ 

n 4>o~ 

Our aim is to prove that A and its adjoint A* defined below are positive in the sense, 

'{au,u) h > 0, yu G V(A), 

(A*U,U) H > 0, \/U G V{A*). 



(U,U) H = (uu + vv + ^&dxdy, \\U\\ H = {(U,U} H } 1/2 . 



(4.4) 
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These properties are needed to apply the Hille-Phillips-Yoshida theorem (see |Yos80tlHP74| ). 
see below. The result for A is now easy thanks to Lemma 4.1 Indeed the following 
calculations are easy, for U smooth in T>{A): 



{AU,U) H = / (u u x + v u y + gcj) x )u + (u v x + v v y + g4> y )v 



+ 



-(u (j) x + v Q (p y + 4>o(u x + Vy))(pdxdy 



(4.5) 



(using integration by parts) 

{—{u +v + —(p ) + g(f>u) 

2 (bn x=0 



dy 



+ 



{^(u* + v* + U*) + g 9 v) 
2 <po 

= h + I 2 + h + h, 

where I\ stands for the boundary term at x = 0, 12 at x = L±, I3 at y = and I4 at y = L2. 
First, the boundary conditions (4.2) imply that I\ = I3 = 0. We then rewrite I2, 1 a as: 

1 



y=L2 



dx 



h 



— [{U 4>) 

2 u 



+ 9 Z (^r-^)v- 2 + v 2 ](L 1 ,y), 



h = — [{V 9) +9 ( 



9 



1 



1 



VQ 



"990 vl 



+ u 2 ] (x, L 2 



which are nonnegative under assumption (4.1). Therefore, we conclude that (AU, U)h > 



for U smooth in T>(A), which is also valid for all U in T>(A) thanks to Lemma 4.1 



We now turn to the definition of the formal adjoint A* of A and its domain T>(A*), in 
the sense of the adjoint of a linear unbounded operator (see |Rud91| ). For that purpose we 
first assume that U 6 *D(A) and U G H are smooth functions, and then compute 



(4.6) 



(AU,U) H = / (u u x + v u y + gcj) x )u + (u v x + v v y + g4>y)v 
Jn 

+ — (u 9x + v 4>y + 90 (u x + v y ))(j)dxdy 
9o 

= Jo + Ji, 

where Jq stands for the integral on Q and J\ for the integral on dVt. For Jq, we have: 



Jn 



(u u x + VQUy + g<f> x )u - (u v x + v v y + g(f> y )v 



- — (uq9x + vo4> y + 9o(u x + Vy))cpdxdy 
90 

= (A*U, U) H , 
where A* is the differential operator defined as follows: 



(4.7) 



_ ~uqu x - VQUy - g<p x 

A*U = I -uqv x - VQVy - g(j) y 

-UQ9X ~ VQ(j)y - (j)0(U X + Vy) J 
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For Ji, taking into account the boundary conditions (4.2), there remains: 

-L 2 



J 







+ 



[u (uu + vv + <t>4>) + g((jm + uc/>)]] (Li,y)dy 
[vq(uu + vv + (jxp) + g(4>v + v(p)~\ {x, L,2)dx. 



According to [Rud91| . V(A*) consists of the U in H such that U h-> (AC/, f/)# is continuous 
on T>{A) for the topology (norm) of H. If {/ is restricted to the class of C°° functions with 
compact support in Q, then J\ vanishes and U ^ Jo can only be continuous if A*U belongs 
to H. If U and A*U both belong to H, the traces of U are well-defined by observing that 



Proposition 2.1 applies to A* as well, and to more general first order linear differential 



operator with constant coefficients. Hence, the calculations in (4.6) are now valid for any 



such U (and U smooth in T>{A)). We now restrict U to the class of C°° functions on which 
belong to T>(A). Then the expressions above of Jo and J\ show that U i-> (AU,U)h can 
only be continuous in U for the topology (norm) of H if the following boundary conditions 
are satisfied: 



(4.f 



0, on r 2 = {x = Li}, 
0, on r 4 = {y = L 2 }. 



Conv ersely, if U,A*U E L (f2) and conditions (4.8) are satisfied, then the calculations 



(4.6) are valid, J\ vanishes, and U h-> (AU, U)h is continuous on T>(A) for the norm of H . 
Therefore we conclude that 



V(A*) = {U eH = L z (ny,A*U E H, and U satisfies (4.8)} 



where we set A*U = A*U, VU E V(A*). 

Let us introduce the density boundary conditions corresponding to (4.8) 

(4.9) U vanishes in a neighborhood of T 2 U T^, 

and define the following space of smooth function with respect to T>(A*): 

V*(0) = {TJ E C°°m 3 , and U satisfies gg}. 

with f = (r 2 u r 4 , r 2 u r 4 , r 2 u r 4 ), we 



Applying Theorem 
obtain 



3.2 



as we did for Lemma 



4.1 



Lemma 4.2. V* (Q) <T\V(A*) is dense inV(A*). 

The proof for the positivity of A* is similar to the proof for A, we thus omit it here. 

5. The mixed case 

In this section, we consider the stationary 2D shallow water equations operator for the 
mixed case. As before, we only consider one case, and thus assume that 



(5.1) 



u < 



v > g4> . 
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Here, we choose the following boundary conditions according to the discussion in Sub- 
section 12.21 



(5.2) 



vqu — uqv + KQ(p = uqu + vqv + g<j) = 0, on Ti 
vqu — uqv — kq4> = 0, on T2 = {x = L\}, 
u = v = 4> = 0, on T 3 = {y = 0}, 



{x = 0}, 



and introduce the corresponding density boundary conditions: 



(5.3) 



vqu — uqv + kq<P, uqu + vqv + g4> vanish in a neighborhood of T\, 
vqu — uqv — Ko4> vanishes in a neighborhood of T2, 
U vanishes in a neighborhood of r 3 . 



We then define the function spaces 



V(A) = {U eH = L 2 (nf,AU G H, and U satisfies Q}, 
V(O) = {U G C°°(U) 3 , and £/ satisfies B}, 



and write At/ = AU, VU G PM). Setting T = (Ti U r 3 ,r 2 U r 3 ,ri U r 3 ), we have that 
P~ 1 U\f = is equivalent to (5.2) and P _1 t7 vanishing in a neighborhood of T is the same 
as (|5.3 ), where P is as in Subsection 2.1 Then we obtain the following result from Theorem 
3.2 with Q = P-^Mi = £i,M 2 = £ 2 ,A = diag(\i,\ 2 ,\ 3 ). 



Lemma 5.1. V(f2) nV(A) is dense in V(A). 



5.1. Positivity of A and its adjoint A*. As we already calculated in (4.5), w e see that 
{AU,U) H = h+h + h + h holds for U G C°°{Tl) 3 . If [7 further satisfies B (i.e. 
17 G V(fi)), we find J 3 = 0, and we rewrite I 4 = f [( v _ JL )2 +5 2 ( _i_ _ ^ + U *J^ ^ 



which is positive by the assumption Vq > 
introduce the new notations £,fj,C,P such that 



It remains to estimate I\ , I2 ■ We first 



/ V -U K \ /A ^ /u\ /v U- U V + K (j)\ 

(5.4) P _1 = v -m -«o , N = \ v \ = \ v ° u ~ u ° v ~ 

\u v g J \(J \(f)J \u u + v v + g<f) J 
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Then we compute I2 as 

\{u,v,4>) 



2/o 




uq g 

u 

q sm 



g 



(5.5) 



(using that fj = v^u 

( uq 
[(|,0,C)P* I no 
5 



M V 



v |](Li,.</) 

— kq4> vanishes in a neighborhood of {x = L\}) 



[( 



u k + gvo 



2 ( u o + v o) K o 



)? + ( 




)C](Li,y) 



> 0. 



Similarly, using that £ = won — + and £ = uou + vow + #0 vanish in a neighborhood 
of {2 = 0}, we compute 2Ii = — [((uqko — 9 v o)/(^{ u o + u o) re o))?7 2 ](0, y) > 0. Therefore, we 
can conclude that (AU, U)h > for C7 G V(fi), which is also true for all £/ G D{A) thanks 
to Lemma I5JJ 



The formal definition of A* can be treated similarly as in Subsection 4.1, we thus omit 
the details. Since we are considering the mix ed c ase, and in order to guarantee that U > 
(AU,U)h is continuous on T>(A) (see Subsec. 4.1), the following boundary conditions need 
to be satisfied: 

' — 

vqu — uqv — kq4> = 0, on ri = {x = 0}, 
(5.6) < vqu — uqv + Ko4> = uqu + vqv + g(f) = 0, on T2 = {x = L\}, 

u = v = 4> = 0, on T4 = {y = L2}. 

We conclude that 

V(A*) = {U G H = L 2 (nf,A*U G H, and U satisfies Q}, 

and write A*U = A*U, W G V{A*). 

We also introduce the corresponding density boundary conditions: 



(5.7) 



Vqu — uqv — KQ(f> vanishes in a neighborhood of Ti, 

vqu — uqv + kq4>, uqu + vqv + g4> vanish in a neighborhood of T2, 

U vanishes in a neighborhood of T4, 



and define the corresponding space of smooth function: 

V*(0) = {U G C(H) 3 , and U satisfies 
Arguing exactly as in Subsection |4.1[ we obtain 



Lemma 5.2. V*(f2) n V(A*) is dense in V{A*), and A* is positive in the sense of (4-4)- 
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6. The fully hyperbolic subcritical case 

In this section, we consider the stationary 2D shallow water equations operator for the 
fully hyperbolic subcritical case. We thus assume that 



3.1) 



ul < g4>o, Vq < g<p , ul + vl > 



Here, we choose the following boundary conditions according to the discussion in Sub- 
section I 



.2) 



vqu — uqv + kqc/) = uqu + vqV + g<p = 0, on Ti = {x = 0}, 
vqu — uqv — Ko4> = 0, on T2 = {x = L\}, 
VqU — UqV — KQ(p = UqU + VqV + g<j) = 0, on T3 = {y = 0}, 
vqu - u v + K cp = 0, on r 4 = {y = L 2 }, 



and introduce the corresponding density boundary conditions: 



(6.3) 



' vqu — uqv + Ko4>, uqu + vqv + g4> vanish in a neighborhood of Ti, 
vqu — uqv — kq4> vanishes in a neighborhood of T2, 
vqu — uqv — kq4>, uqu + vqv + g(f) vanish in a neighborhood of T3, 

k vqu — uqv + kq4> vanishes in a neighborhood of T4. 



We then define the function spaces 



V(A) = {U G H = L 2 (n) 3 ,AU £ H, and U satisfies 
V(O) = {U G C°°(n) 3 , and U satisfies B}, 



and write AU = AU, VU E Setting T = (Ti U r 4 ,r 2 U r 3 ,ri U r 3 ), we have that 

P~ l U\f, = is equivalent to (6.2) and P~ 1 U vanishing in a neighborhood of T is the same 
as (|6.3 ), where P is as in Subsection 2.1 Then we obtain the following result from Theorem 
3.2 with Q = P- 1 , Mi = S 1} M 2 =S 2 ,A = diag(X u A 2 , A 3 ). 



Lemma 6.1. V(£l)nV(A) is dense in V(A). 



6.1. Positivity of A and its adjoint A*. As indicated in Subsection 5.1, we have that 
(AU, U) H = h+h + h + h holds for U G C°°(n) 3 . Suppose that £7 belongs to V(O). In 
order to estimate I±, I 2 , 13, h, we still use the notations £,f},£,P introduced in Subsection 



5.1 (see (|5.4|)), and I\, I 2 are both nonnegative since the estimates for I\,I 2 are exactly the 
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same as in Subsection |5.1[ Proceeding exactly as for I2, we compute I4 as: 



2L 



[(u,v,(f>) 



[foMP* I 

g 




v |](ar } £, 2 ) 



3.4) 



(using that £ 



vqu — uqv + kq4> vanishes in a neighborhood of {y = L 2 }) 







'0 N 



«o 

9 \P\fj \}(x,L 2 ) 

q „ gvo 



g 



4>o 



v k + gu , 2 



"0 



;)C 2 ](^2 







> 0. 



Similarly, using that 77 = vo"U — uqv — kq and ( = uqu + uo^ + 50 vanish in a neighborhood 
of {y = 0}, we compute 2/3 = — [((fo^o ~~ g u o) / (2{ u o + ^o^o))?? 2 ]^) 0) — 0- Therefore, we 
conclude that (AU,U) H > for U G V(Q); this implies that (AU,U) H > also holds for 
all f7 € f(-A) by virtue of Lemma |6.1| Hence A is positive. 



Taking similar arguments in Subsection 4.1 and Subsection 5.1 we will obtain the same 
results for A* the adjoint operator of A. We thus only state them below without the proof. 

In the fully hyperbolic subcritical case, we first introduce the following boundary condi- 
tions , 



3.5) 



vqu — uqv — kq4> = 0, on r± = {x = 0}, 
vqu — uqv + KQCp = uqu + vqv + g<t> = 0, on T2 = {x = L\}, 
vqu - uqv + Ko<fi = 0, on r 3 = {y = 0}, 
{v u - u v - K (j) = u u + v v + g(f> = 0, on r 4 = {y = L 2 }; 

and the corresponding density boundary conditions: 

vqu — Uqv — kq4> vanishes in a neighborhood of F\, 
vqu — uqv + kq(J), uqu + vqv + gcj) vanish in a neighborhood of T 2 , 
vqu — uqv + HQ(f) vanishes in a neighborhood of T3, 
^vqu — Uqv — kq(J), uqu + vqv + g<j) vanish in a neighborhood of 1^4. 

We then define the function spaces: 

V{A*) = {U G H = L 2 (n) 3 ,A*U G H, and U satisfies (JO}}, 

V*(fi) = {U £ C{Tl) 3 , and U satisfies Q}, 

and write A*U = A*U, V£7 G V(A*). Then we have the following results. 



3.6) 



Lemma 6.2. V*(0) n V{A*) is dense in V{A*), and A* is positive in the sense of (4-4)- 
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7. A PRELIMINARY REGULARITY THEOREM 

We aim in Section [8] to study the mixed subcritical case which corresponds to 
u l + v l < 94>0, implying u 2 < g<j) , Vq < gfo, 

and which embodies some elliptic operators. In this section we introduce some necessary 
preliminary results, namely we prove a general regularity theorem regarding functions de- 
fined on the domain f2 = (0, L\) x (0, L2). This theorem will be very useful in verifying the 
validity of the integration by parts in Section |8j Following the notations in Section [3j we 
introduce the following boundary conditions: 



(7.1) 



on r, 
on r c 



where T c is the complement of T with respect to the boundary d£l. Here we choose T 
Ti U T3, which is the case considered in Section [8} We define 

V = {9 = (61,62)* G H 1 ^) 2 I O satisfies (H) }. 

Note that ||VS|| i2 (Q)2 is a natural norm on V thanks to the Poincare inequality. 
We assume that ai, «2, A., P2 are four real constants satisfying 

(7.2) avPi - axfa ± 0, 
and we define the operator T such that 

(7.3) re = T 1 o x . + r 2 ?/ , 

where 

rp foil Pl\ ^ _ (0L 2 P2 



/3i —aij ' \^2 —OL2 J 

Notice that T is symmetric and elliptic, and actually we have the following statement. 

Proposition 7.1. ||'7 r " ©|| i2 (o) 2 * s a norm on V equivalent to the H^-norm. 

Proof. Let us recall a basic fact from linear algebra. We endow R 2 with its usual dot product 
and induced norm: 

x ■ y = xiyi + x 2 y 2 , \x\ 2 = \Jx\ + x?,, 

where x = (xi, x?), y = (2/1,2/2) G R 2 - Let T be a linear transformation from M? to IR 2 with 
To as its matrix representation, where 

(oil a 2 

T °-Ui h 



Condition (7.2) shows that To is non-singular, which implies that T is an isomorphism. 
Hence we have that 

(7.4) — \x\ 2 < \Tx\2 < c 2 \x\2, 

Cl 

where C2 (resp. ci) is the spectral norm of To (resp. Tq 1 ), i.e. the square root of the largest 
eigenvalue of the positive-definite matrix TqTq (resp. T _ *T _1 ). 
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We now compute 



Ire! 



L 2 (n) 2 



(7.5) 



l (0*2? + ©^(T^ + T 2 Q y )dxdy 
Jn 

jf (a? + /3?)(C + O + («! + + 

+ 2(aia 2 + PifoWixOiy + 2x 6 2y ) 

+ 2(a 2 Pi - ai(3 2 )(9 2x 9i y - Q lx Q 2y )dxdy. 



ly) 



To dispense with the last term in the integral of ( |7.5[ ), we use a result from [Gri85j (see 
Lemma 4.3.1.3), which implies: 



Lemma 7.1. The identity 



>2xV\y 



dxdy 



hx02ydxdy 



holds for all © = (0i,0 2 )* G i/ 1 ^) 2 satisfying (7.1) (i.e. OeV). 



Thanks to Lemma 7.1 (7.5) gives 

f (aie lx + c^) 2 + {Pi9i x + fo6i y f 
Jn 



TO 



L 2 (Q) 2 



(7.6) 



n 



+ (ai0 2 x + o^) 2 + (/?i02x + (3 2 e 2y ) 2 dxdy 
TV0i|! + iTV^alldxdy. 



With the help of (7.4), (7.6) immediately implies that 
1 

ci 



Hence Proposition |7.1| follows. 

The regularity result is the following existence theorem. 



□ 



Theorem 7.1. For every given F = (/i, f 2 ) t G L 2 (f2) 2 , the problem 7"@ = F has a unique 
solution € V. 

Proof. Let /(©) = (F,TQ) and a(@,@) = (TQ,T@) L 2 (n) 2 ^ then it is not hard to check 
that the linear form I is continuous on V, and the bilinear form a is continuous on V x V. 
Observing that the form o is coercive on 1/ thanks to Proposition 7.1 we obtain that there 
exists a unique G V such that 

(7.7) a(O,0) = /(©), 

for all G V thanks to the Lax-Milgram Theorem. 

In order to interpret (7.7) in the distribution sense, we need the following lemma. 

Lemma 7.2. For every given ^ = (ipi,^) 1 G V(Q) 2 , the problem T0 = \l/ has a unique 
solution Q £ V . 
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7.1 



The proof of Lemma |7,2| will be given below, and we continue with the proof of Theo rem 
For any \& G V(Q) 2 , we find G V such that TO = \& by virtue of Lemma 7.2, and 

F 

□ 



then (7.7) shows that TO = F in the distribution sense in Q, and we conclude that TO 
in L 2 (?Tp since V(n) 2 is dense in L 2 (VL) 2 . 

Proof of Lemma \ 7. S\ Without loss of generality, we first assume that, 
(pl) a 2 fa - a x fa = 1. 



Let ^ = (V ; i)'02) < £ T>(£1) 2 : we look for G V satisfying TO = \£. In order to find such 
O G V, we introduce a new coordinate system (a/, y') such that 



ft 

«2 



"A 

-ai 



which is equivalent to 



-u\ ft 
-a 2 ft 



fax - fay 
a 2 x — a±y 



-a\x' + fay' 
-a 2 x' + fay' 



We denote by the image of Ti by this transformation for alii G {1, 2, 3, 4}, and denote by 
f2', r', 9', *&' and the gradient V' the transforms of f2, T, 9, and the gradient V respectively. 
Now, direct computation gives 



(7.8) 



V9 



fa 
-fa 



a 2 



V'9'. 



In the new coordinate system (x',y'), the boundary conditions (7.1) read 

<E3) 



?i = on r', 
?i=0on r ,c , 



where T' = T[ U r' 3 . We also denote by V the function space corresponding to V: 

V' = {©' = (9' 1 ,9' 2 ) t G H 1 ^') 2 | ©' = (9' 1 ,9' 2 ) t satisfies (TtJ }. 
In the new coordinate system (x' , y'), the operator T reads 



TO' 



fa fa 



V0i + 



(7.9) 



(using (7.8[)) 

ft ft 



ft 
-ft 



ft 
-ai 



«2 

— ai 



ft 

-«2 



V0 2 



ft 

-ai 



ft 

-«2 



ft 
"ft 



0:2 
— Ql 



(using assumption (7.2 )) 



-1 



1 
ol 



V'9[ + 



1 

1 



V9' 2 



'lx 



, + 



2y' 

2x' 
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and of course, we have VP' 6 T>{0!). We are seeking ©' G V satisfying T@ = that is 
(7.10) 



+J'2y> = ^'l, 
2x> ~ 



Differentiating (7.10)i with respect to x' and (7.10)2 with respect to y' ', and subtracting 
these two equations, we find the elliptic equation 

(7.11) tie' l = -i>' lx , + i>' 2y „ 

where A' denotes the Laplace operator in the new coordinate system (x',y'). We associate 
to equation (7.11) the boundary conditions 

(7.12) 9[ = on r, 



which is already contained in (7.1 ]). For the boundary r /c , suitable boundary conditions 
can be obtained as follows. We first denote by v[ the unit normal vector to 1^, and r- the 
unit tangent vector on T^, for all i £ {1,2,3,4}. On T' 2 = {— a\x' + Piy' = Li}, we have 
d9' 2 /dr 2 = since 9' 2 = on T 2 . Noticing that (Pi, a\) is parallel to t 2 , we find 



(7.13) 



= (f3x,a 1 ) 



y 



on r' 2 . 



Multiplying (7.10)i by a>i, (7.10) by Pi, and adding these two equations, we find 
- ai6' lx , + PiO'iyi + ai9' 2y , + P\9' 2x , = aiip^ + /3iV4> 



which, by using (7.13) and noticing that ^' = (ip'i,^) vanishes on T 2 , implies that 



(7.14) 



-ax, Pi 



y 

'lx' 

y 



e' lx , + Pi0' ly , = o, on r' 2 . 



Observing that v' 2 is parallel to (— a\, Pi), we then associate to (7.11 ) the following boundary 



condition on T' 2 : 



(7.15) 



d9[ 



0, on r' 2 , 



which is equivalent to (7.14). On T' 4 = {— a^x' + P^y 1 = L2}, we have d9 2 /dT' 4 = since 
9' 2 = on r' 4 . Noticing that (/^c^) is parallel to T4, and u' A is parallel to {—0.2^2)-, similar 
computations show that we need to associate to ( |7.11[ ) the following boundary condition on 



V ■ 

1 4- 



(7.16) 



d9[ 
di7 A 



0, on r^. 



The existence and uniqueness of a solution 9[ £ H 1 ^') of (7.11)-(7.12) and (7.15)- 
( [77l6| ) follows from Lemma 4.4.3.1 of |Gri85| with Sl = Q,',& = T^yV ^T^,Pj = 0,/ = 

-1p lx , +lp 2 y>. 
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Following the arguments for 9{, we find the following equation and boundary conditions 
for 9 2 : 



'A'O! 



(7.17) 



2 - fly' +^2x>, 

2 = o, on r c = r 2 u r' 4 , 



du_[ 



o, on ri, 
o, on r' 3 . 



We thus also have a unique solution 0' 2 G H 1 ^') of (7.17) thanks to Lemma 4.4.3.1 of 
[Gri85] again. 

In conclusion, in the new coordinate system (x',y r ), we find a unique solution O that 
belongs to V' , and solves the problem TO = Transforming back to the original coor 
dinate system (x,y), we obtain O 6 V satisfying TO = Hence, the proof of Lemma 7.^ 
is complete. 

□ 



Remark 7.1. Theorem 7J_ is still true, if we choose any other boundary T except and 
U^ =1 Fj. The proof is exactly the same. 

We will now develop a similar result for T* the adjoint operator of T. The adjoint T* is 
formally calculated as usual. For smooth functions O and O, we compute 



TO • edxdy 



©ViOz + e t T 2 Q y dxdy 



(using integration by parts) 



@ t T l Q x + Q^Qydxdy 



+ 



(7.18) 



O TO 



x=Li 



x=0 



dy + / ©TO 



V=L 2 
y=0 



dx 



(using the boundary conditions ( |7.1| )) 
T*0 • Qdxdy 

L 2 

(a x h + hh)0i 

x=L 



+ 



{a 2 h + p 2 e 2 )6 y 



y=L 2 



{Mi ~ »i9 2 )9 2 
(M - a 2 9 2 )9 2 



x=0 



dy 



dx, 



where T*0 = —T\Q X — T 2 @ y = — T@. Hence, we introduce the following boundary condi- 
tions corresponding to the adjoint operator T* (dropping the bars): 



(7.19) 



r Mi-ai9 2 


= 0, 


on 




= {x 


= 0}, 


{ 0101 + 0102 


= 0, 


on 


r 2 


= {x 


= L 1 } 


(3 2 9l — «2#2 


= 0, 


on 


r 3 


= {y 


= 0}, 


a 2 9i+p 2 9 2 


= 0, 


on 


r 4 


= {y 


= L 2 } 
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and the function space 



V = {Q = (0i, 2 ) G H L (ny I 9 satisfies fl7l9[ ) }. 

The results for T* associated with V are the same as for T associated with V. We state 
them as follows. 

Proposition 7.2. ||7^*0| l^a^a ^ s a norm on V equivalent to the H l -norm. 



The proof of Proposition |7.2| is the same as the p roof o f Proposition |7. 1| since Lemma 7.1 
also holds for all 6 = (0i,0 2 )* G H^fl) 2 satisfying gig} (i.e. 6 G F). 

Theorem 7.2. For every given F = (/i, / 2 )* G _L 2 (Q) 2 , i/ie problem T*@ = F has a unique 
solution G V. 



The proof of Theorem 7.2 is also the same as the proof of Theorem 7.1 using the following 



lemma which replaces Lemma 7.2 



Lemma 7.3. For every given ^ = (ipi,ip2) t G V{Q) 2 , the problem T*@ = \P possesses a 
(possibly non-unique) solution G V. 

Proof. To show the result, we introduce a suitable change of variables (different from the 
previous one) which makes the boundary conditions (7.19) simpler. We set 



"101 + /?10_2 

P\Qi - a±02 



ai Pi 
Pi -«i 



Ti6, 



and find by direct computation 

/3 2 0~i - a 2 2 = Ml£i + O201 + /3202 = M26 - Mi6> 
where the constants fi± and /x 2 are given by 

OC1P2 - OL2P1 aifii + a 2 /3 2 



Ml 



M2 



Thus the boundary conditions (7.19) are equivalent to 

'6 = 0, onr 1 = {x = 0}, 
£1 = 0, on r 2 = {x = Li}, 
Ml£l + M26 = 0, on r 3 = {y = 0}, 
.M26 - M16 = 0, on r 4 = {y = L 2 }. 



(7.19') 



With the new variables, we define the operator V such that VE = 7~*0, that is 



VE = T*@ = -TiQ x - T 2 



(7.20) 



£1 

1 M2 

/xi 



y — ^ 

M2 "Ml 

Mi M2 



£1 
£2, 



V£i 



-Mi 

M2 



V£ 2 . 
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The goal now becomes the following: given ^ = (ipi,^) 1 G ~D(Q) 2 , we look for E G V such 
that = ^, where the function space V is defined as 



F = {E = (6,6)* G tf 1 ^) 2 I 3 satisfies (7.19') }. 



Noticing that /ii is nonzero because of assumption (7.2), we can assume without loss 
of generality that \x\ = 1. Following the proof of Lemma 7.2, we then introduce a new 
coordinate system (x',y') such that 



x 



1 



1 

M2 



We denote by V the gradient in the new coordinate system (x',y') In the new coordinate 
system (x',y'), the gradient V and the operator V read 



(7.21) 



-M2 
1 



v'6 



(7.22) 



/i 2 
1 

1 



"A*2 

1 



-1 

/'2 



V6 



-l 

/<2 



"M2 
1 



V'6 



(7.19 1 ) 



and the boundary conditions (7.19') read 

6 = 0, onTi, 
6 = 0, on r 2 , 
6 + ^26 = 0, on r' 3 , 
A«26 -6=0, on r 4 . 
Finally we define the function space 

V i = {E = (6,6)* G tfW) 2 I 3 satisfies ( |7.19f| ) }. 
We are now seeking HgV satisfying = that is 

-6a' +6y = i>i, 
-6y - 6x' = ip2- 



(7.23) 



We now have two cases to consider. First, if /i 2 = 0, then the boundary conditions (7.19') 



are simpler, and similar to the boundary conditions (7.1 ). Hence, we obtain a unique 

solution S G V that solves the problem "PH = ^ by following the same argument as for 
Lemma 7.2. If /i 2 ^ 0, then without loss of generality, we can assume that = 1, and the 



arguments presented below will be similar to the proof of Lemma 7.2 



Differentiating (7.23)i with respect to y' and (7.23)2 with respect to x' , and subtracting 

A'6 = ipltf - fax' • 



these two equations, we find 
(7.24) 
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We first associate with equation (7.24) the boundary condition 



(7.25) 



6 = o, on ri, 



which is already contained in (7.19'). For the other boundaries, suitable boundary con- 



ditions can be obtained as follows. As in Lemma 7.2, we denote by v\ the unit normal 



vector to 1^, and t[ the unit tangent vector on 1^, for all i € {1,2,3,4}. On T' 2 , we have 
d^i/dr^ = since £i = on T 2 . Noticing that (0, 1) is parallel to t 2 , we find that = 0, 
which, together with (7.23)2, implies that ^2x' = — Ciy' = by noticing that ip2 vanishes 



on T' 2 . Observing that v 2 is parallel to (1,0), we then associate to (7.24) the following 
boundary condition on T 2 : 

(7.26) S = 0. 



dv' 2 



on T 2 . 



Similar arguments show that we need to associate to (7.24) the following condition on the 
other two boundaries: 



(7.27) 



du' 
du' 



+ 



Or' 
Or' 



0, on H 



3' 



o, on r 4 . 



The existence of a (possibly non-unique) solution £2 £ H l (Cl') of ( 7.24 )-( 7.27) follows from 
Lemm a 4.4.4.2 of |Gri85| with Q = ft', @ = T[, jV = r^Ul^U! ^, / = ipi y '-ip2x>- Integrating 
(7.24)2 with respect to y 1 and using the boundary condition (7.19 '(3, we obtain a solution 



£1 which belongs to H l (Q!). 

In conclusion, in the new coordinate system (x',y'), we find a (possibly non- unique) 
solution H that belongs to V , and solves the problem VE = ^> . Transforming back to the 
original system (x,y) and the original variables O, we obtain 9 6V satisfying 7~*0 = \P. 
Therefore, the proof of Lemma 7.3 is complete. □ 



8. The mixed subcritical case 

In this section, we consider the stationary 2D shallow water equations operator for the 
last case which includes both hyperbolic and elliptic modes. We assume that 

(8.1) ul+vl < gfo, 
which implies that 

(8.2) ul < #0, vl < gfo. 

Here, we choose the following boundary conditions according to the discussion in Subsection 



vou — uqv = uqu + vov + gcj) = 0, on ri U T3 = {x = 0} U {y = 0}, 
= 0, on T 2 U T 4 = {x = Li} U{y = L 2 }. 



(8.3) 

We then define the function space 



V(A) = {U £H = L 2 (nf,AU £ H, and U satisfies gjj}, 
and set AU = AU, VU G T>(A). 
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8.1. Positivity of A and its adjoint A*. We use the notations £, rj, £, P, So introduced 
in Subsection 2.4 and set H = (£, 77, so that [7 = PH. We then write 



Ai = P'SoSiP 



A 2 = P t S £2P 



and for U £ D(A), we compute: 




ul + vl 



(AU, U) H 



n 



(u u x + v u y + gcf> x )u + (u v x + + g<t> y )v 



9 



i.4) 



+ -^-{uq^ x + uo^y + (/>o(mx + v y ))(f>dxdy 

(S £iU x + S £ 2 U y ,U) H 
(ptSoSiPZ* + P t S £ 2 PZ y , S) H 



+ 



u 



gyp 

Kl 

-n 



gup 

Kl 



guo 
-VQ 



dxdy 



ul + vl J n 



COoCx + v ( y )dxdy. 



We now need to apply integration by parts for the two terms Jo, J\ in the right-hand side 
of (8.4). According to Theorem 7.1, we obtain that £,77 actually belong to H l {^l) since £,77 
satisfy the elliptic equations (2.17). Hence, integration by parts are valid for Jo, and we 
find 



Jo 



2(«§ + vl) Jo 
+ 



mi 



gyp 

Kl 

-u 



x=L± 



1 



2K + «o) 



Vo 
gup 

Kl 



x=0 

gup 

Kl 

■Vo 



dy 



y=L 2 



dx. 



y=0 



For Ji, with density Theorem 3.1 applied with 9 = (, A = uo/^o,r = Ti U T3, we see that 
the integration by parts are also valid for Ji, and we find 



Ji 



2{u 2 + vl) 



uoC 



x=L\ 



x=0 



dy + 



voC 



y=L 2 

y=0 



dx 



Combining the results for Jq and Ji, we arrive at (written in a compact form): 



i.5) 



(AU,U) H 



L 2 



x=Lx 1 

dy + 7y 
x=o 2 



~<A 2 ~ 



y=L 2 
y=0 



dx 



h+h + h + h, 



30 



A. HUANG AND R. TEMAM 



where Ii stands for the boundary term at Tj for i = 1, 2, 3, 4. Using the boundary conditions 
(8.3), we now estimate Jj, (i = 1, 2, 3, 4) as follows: 

2{u 2 + v 2 )h = -(u f + 2^£r l -u ri 2 + u ( 2 )(0,y) = (u oV 2 )(0,y) > 0, 

2{u 2 + vl)I 2 = (not 2 + 2^7? - uoif + 7z C 2 )(Li, y) = K£ 2 + n C 2 )(0, y) > 0, 

>3 u o, 



2(^ + Wo 2 )/ 3 = -(^ 2 



2^7, - W + «oC)(x, 0) = (v oV z )(x, 0) > 0, 

Kl 

,5^0 2 , /-2w r ^ / t 2 , >2\ 



2(t^ + ^)J 4 = (yot ~ 2— ^ - W + voC)(x, L 2 ) = (vq? + v ( 2 )(x, L 2 ) > 0, 

Kl 

which are all nonnegative. Therefore, we conclude that (AU,U)h > for all U E T>(A). 
Note that here, because the part of the operator A which is elliptic gives some additional 
regularity, we did not need to use an approximation argument to show that (AU,U)h > 
0, \/U G D(A), unlike in the previous cases. 

We now turn to the adjoint A*. In the mixed subcritical case, the formal definition of 
A* can also be treated similarly as i n Su bsection 4.1, we thus omit the details. The adjoint 
differential operator A* is given by (4.7) again. In order to guarantee that U h-> (AU,U)h 
is continuous on T>(A), the following boundary conditions in the (£, rj, £) variables need to 
be satisfied (see (7.19) since the equations for £,77 are elliptic): 

(8.6) 



f ^-u f} = 0, onTi = {x = 0}, 
C = ff 77 = 0, on r 2 = {x = Li}, 

- 9 ^-v fj = 0, onr 3 = {y = 0}, 

0, on r 4 = {y = L 2 }. 



Transforming back to the (u, v, 4>) variables, the boundary conditions read 

gv^u — qvqUqv — uon 2 4> = 0, on Fi = {x = 0}, 
uqvqu — UqV + gvo4> = uqu + vqv + g<j) = 0, on T 2 = {x = L{\, 



(8.7) 



gu v u - gu\v + vqk\4> = 0, on r 3 = {y = 0}, 



vlu - vqUqv - gu (p = u Q u + v v + 



0, on r 4 = {y = L 2 }, 



which are more complicated than the boundary conditions (8.3), since some elliptic modes 
are hidden in the operator A. We conclude that 

V(A*) = {U G H = L 2 (nf,A*U G H, and U satisfies Q}, 

and set A*U = A*U, \/U G V(A*). 

The proof of the positivity of A* is similar to that of A, we thus omit it here. 



9. The initial and boundary value problem for the full system 

In this section we aim to combine the results of the previous sections and to investigate the 
well-posedness for Eqs. (1.1) associated with the suitable initial and boundary conditions 
that we have introduced. Here again, we consider the case of homogeneous boundary 
conditions (see below a remark about the non- homogeneous boundary conditions). 
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As we have already seen in previous sections, we have five cases to consider depending 
on the relations between the constants Uo,Vo,(j)o,g. In the following, we list these relations 
and their corresponding boundary conditions: 

The supercritical case. The assumption is 



(9.1) u > g<p , v Q > gcpo, 

and the corresponding boundary conditions are 



u = v 

U = V 



0, on T 1 = {x = 0}, 
0, on r 3 = {y = 0}; 



The mixed case I. The assumption is 

(9.2) ul < g(f) , Vq > g4> , 
and the corresponding boundary conditions are 

vqu — uqv + KQtfi = uqu + vqv + g<p = 0, on Ti = {x = 0}, 
( |9.2| ) < vqu — uqv — kq(J) = 0, on T2 = {x = L\}, 

u = v = (j) = 0, on T3 = {y = 0}; 

The mixed case II. The assumption is 

(9.3) ul > g<p , Vq < gfa, 
and the corresponding boundary conditions are 

v u - u v - kq4> = uqu + vqv + gcf> = 0, on T 3 = {y = 0}, 
( |9.3| ) < v u - u v + K (p = 0, on T 4 = {y = L 2 }, 

u = v = 4> = 0, on Ti = {x = 0}; 

The fully hyperbolic subcritical case. The assumption is 



(9.4) 



u < 



, Vq < g<f) , u 2 + vl> 



and the corresponding boundary conditions are 

vqu — uqv + KQ(p = uqu + vqv + g<j) = 0, on T\ = {x = 0}, 
vqu — uqv — KQ(p = 0, on T2 = {x = Li}, 
vou — uqv — Kocp = uqu + vqv + g<f> = 0, on T3 = {y = 0}, 
{v u - u v + K <j> = 0, on r 4 = {y = L 2 }; 

The mixed subcritical case. The assumption is 

(9.5) 

and the corresponding boundary conditions are 



Uq < 94>o, < gfa, v% + Vq < 



dpi) 



vqu — uqv = uqu + vqv + gcj) = 0, on Ti U T3 = {x = 0} U {y = 0}, 
= 0, on T 2 U T 4 = {x = Li} U {y = L 2 }. 
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If these constants Uq, vq, (f>o, g satisfy assumption (9.i), then we define the domain of the 
unbounded operator A: 

V(A) = {U G H = L 2 (nf, AU G H, and U satisfies (9.i)}, 

where i G {1, 2, 3, 4, 5}, indicating the five different cases. 

We set BU = (—fv,fu,0Y, where / is the Coriolis parameter, and it is easy to see 
that B is a linear continuous self-adjoint operator on H. We also set Aq = A + B, with 
V(Ao) = T>(A), and the adjoint A% = A* + B*, with V(A* ) = V(A*). Then we have 

Theorem 9.1. The operator —Aq is the infinitesimal generator of a contraction semigroup 
on H. 

Proof. According to [Yos80, HP74], it suffices to show that 

(i) Aq and Aq are both closed operators, and their domains T>(Aq) and T>(Aq) are both 
dense in H. 

(ii) Aq and Aq are both positive. 

Noticing that B is a linear continuous self-adjoint operator on H and that {BU, U) = for 
all U G H, we find that proving (i) and (ii) is equivalent to proving the following: 

(f) A and A* are both closed operators, and their domains T>(A) and D(A*) are both 

dense in H. 
(ii') A and A* are both positive. 

We already showed that A and A* are both positive in the previous subsections, we thus 
only need to prove (f). We establish the result for A, and the proof for A* would be similar. 

In all cases, observing that P(f2) 3 is included in T)(A) and dense in H = L 2 (S1) 3 , we 
thus obtain V{A) is dense in H. To show that A is closed, consider a sequence {U n } ne n C 
V(A) for which lim n ^oo U n = U in L 2 (0) 3 and lim^oo AU n = U in L 2 (ft) 3 . By the L 2 - 
convergence of {C/ n } ng N, we first find that AU n converges to AU in the sense of distributions, 
which implies that AU equals U in the sense of distributions. Since U belongs to L 2 (S1) 3 , 
we obtain that AU belongs to L 2 (f2) 3 , too. Thus, we have U,AU G L 2 (fi) 3 , which shows 



that the traces of U are well-defined thanks to Proposition 2.1 By Proposition 2.1 again, 
the traces of U n converge weakly to the traces of U in the appropriate space H~ l , so that 
U satisfies the boundary conditions (9.i'), where i G {1,2,3,4,5}. Therefore, we conclude 
that U G D(A). Hence A is closed, and the proof is complete. □ 

We now consider the whole system of 2D linearized Shallow Water Equations, namely 



(1.1) and introduce the initial and boundary conditions. We consider the homogeneous 
boundary conditions and will briefly explain how to treat the case of non-homogeneous 
boundary conditions. 

As we discussed before, the boundary conditions are (9.i'),i G {1, 2, 3, 4, 5} depending on 
the case we are considering, and all these boundary conditions are taken into account in 
the domain T>(Aq) of Aq. Finally if we add the initial conditions: 



(9.6) 



U(0) = («(0), 17(0)^(0)) = U = (uq,vq,cPo), 
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then the initial and boundary value problem consisting of Eqs. ( 1.1 ), (9.i'),i G {1, 2, 3, 4, 5} 
and (9.6) is equivalent to the abstract initial value problem 

(97) f% + AoU = F, 



Note that F = (F u , F v , F^) which does not appear in (1.1) is added here for mathematical 
generality and to study the case of non-homogeneous boundary conditions. Thanks to 
Theorem 9.1 this problem is now solved by the Hille-Yoshida theorem and we have that: 

Theorem 9.2. Let H,Aq and T>(Aq) be defined as before. Then the initial value problem 
(9.7) is well-posed. That is, 

i) for every Uq G H, and F G L 1 ^, T; H), the problem (9.1) admits a unique weak solution 
U G C([0,T];H) satisfying 

ft 



U(t) = S(t)U + S{t- s)F(s)ds, Vt G [0, T], 
Jo 

where (S(t))t>o is the contraction semigroup generated by the operator —Aq; 
ii) for every U ~e V(A ), and F G L^TjH), with F' = dF/dt G L l (0,T;H), the 
problem (9.1) has a unique strong solution U such that 

UeC([0,T];H)nL°°(0,T;V(A )), ^ G L°° (0, T; H) . 

Remark 9.1 (Non-homogeneous boundary conditions). In the case of non-homogeneous 
boundary conditions, that is if we want to solve (1.1), with (9.H), i G {1,2,3,4,5}, in which 
the boundary conditions are now non-homogeneous, and with initial condition (ff.6| ) ; we 
classically proceed by translation. We assume that the boundary data are inferred from a 
function $ which is defined on Q x [0,T], and we set V = U — The problem is now to 
check the hypotheses of Theorem\9. S\ for the data corresponding to V and to apply Theorem 



9.2. We omit the details; see [RTT08bj for a similar situation. 



Remark 9.2. The linearized 3D inviscid Primitive Equations of the atmosphere and the 
oceans can be written: 



(9i 



u t + U Q u x + V u y - fv + <j) x 
vt + U v x + V v y + fu + 4> y 
T t + U T x + V T y + N 2 fw 



U X +Vy+ W z 

Yz — t ' 



■V 

0, 



o, 

0, 



where (u, v, w) is the velocity of the water, (u, v) the horizontal velocity, T the temperature, 
4> a multiple of the pressure; g the gravitational acceleration, and N 2 denotes the Brunt- 
Vdisdld (buoyancy) frequency satisfying 

_g_dp 
Po dz' 

and Uq > 0, Vq > 0, po > and Tq > are reference average values of the density and the 
temperature. 
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We refer to |RTT08b] for a systematic discussion in the case when Vq = 0. Using the 



normal modes expansion for (9.8) (see details in [RTT08b] ), we obtain the equations for 
the non-zero modes (u n , v n , <f) n , A n ) which read (the indices n are dropped for the sake of 
simplicity): 

u t + U u x + V u y - fv - jip x = 0, 
(9.9) I v t + U v x + V v y + fu- \i> y = 0, 

^ t + U 1p X + V 1py ~ ^(U X + Vy) = 0, 

where ip = <fi z and X is positive. Observing t hat i9.9\ ) has the same form as the shallow 
water equation (1.1) if we replace <f> and g in (1.1) by —ip and A -1 , we thus can obtain a 
well-posedness result (see Theorem 9.2) for \9.ty . The equations for the zero mode are 

u t + U u x + V u y - fv + \<j} x = 0, 



(9.10) 



Vt + U Q V X + V Vy + fU + J(j)y = 0, 



u x + v % 



0. 



Following the same argument as for the zero mode in [RTT08b| and [CST10j ; we can also 
obtain the well-posedness result for (9.10). Therefore, we can obtain a well-posedness result 
for the whole system (9.8) with suitable boundary conditions and initial conditions as in 
Section 4 of [RTT08b] . The details will appear elsewhere. 

Acknowledgments 

This work was partially supported by the National Science Foundation under the grants 
NSF DMS-0906440 and DMS-1206438, and by the Research Fund of Indiana University. 



References 

[AB03] E. Audusse and M.-O. Bristeau, Transport of pollutant in shallow water: A two time steps kinetic 
method, M2AN Math. Model. Numer. Anal. 37 (2003), no. 2, 389-416. 

[ABBKP] E. Audusse, F. Bouchut, M.-O. Bristeau, R. Klein, and B. Perthame, A fast and stable well- 
balanced scheme with hydrostatic reconstruction for shallow water flows, SIAM Journal of Scientic 
Computing 25 (2004), no. 6, 2050-2065. 

[ABPS] E. Audusse, M.O. Bristeau, B. Perthame, and J. Sainte-Marie, A multilayer saint-venant system 
with mass exchanges for shallow water flows, derivation and numerical validation, M2AN Math. 
Model. Numer. Anal. 45 (1) (2011), 169-200. 

[BS07] S. Benzoni-Gavage and D. Serre, Multi-dimensional Hyperbolic Partial Differential Equations, 
Oxford University Press, 2007. 

[BPSTT] A. Bousquet, M. Petcu, M.-C. Shiue, R. Temam, and J. Tribbia, Boundary conditions for limitied 
area models, To appear. 

[CST10] Q. Chen, M.-C. Shiue, and R. Temam, The barotropic mode for the primitive equations, Speical 
issue in momory of David Gottlieb, Journal of Scientific Computing 45 (2010), 167-199. 

[Eva98] L. C. Evans, Partial Differential Equations, vol. 19, Amer. Math. Soc, Providence, RI, 1998. 

[Gri85] P. Grisvard, Elliptic Problems in Nonsmooth Domains, Monographs and Studies in Mathematics, 
Pitman, Boston, 1985. 

[Hor65] L. Hormander, L? estimates and existence theorems for the 3 operator, Acta Math. 113 (1965), 
89-152. 

[HPT11] A. Huang, M. Petcu, and R. Temam, The one- dimensional supercritical shallow-water equations 
with topography, Annals of the University of Bucharest (Mathematical Series) 2 (LX) (2011), 
63-82. 



THE LINEARIZED 2D INVISCID SHALLOW WATER EQUATIONS 



35 



[HP74] E. Hillc and R.S. Phillips, Functional Analysis and Semi-Groups, American Mathematical Society, 
Providence, HI, 1974 (Third printing of the revised edition of 1957, AMS Colloquium Publications, 
vol. XXXI). 

[Kre70] H.-O. Kreiss, Initial boundary value problems for hyperbolic systems, Comm. Pure Appl. Math 23 
(1970), 277-298. 

[KT80] K. Kojima and M. Taniguchi, Mixed problem for hyperbolic equations in a domain with a corner, 

Funkcialaj Ekvacioj 23 (1980), 171-195. 
[Lop70] Ya. B. Lopatinskii, The mixed Cauchy-Dirichlet type problem for equations of hyperbolic type, 

Dopovfdf Akad. Nauk Ukrai"n. RSR Ser. A 668 (1970), 592-594. 
[OS78] J. Oliger and A. Sundstrom, Theoretical and practical aspects of some initial-boundary value 

problems in fluid dynamics, SIAM J. Appl. Math. 35 (3) (1978), 419-446. 
[Ore95] P. Orenga, Un th'eoreme d'existence de solutions d'un probl'eme de shallow water, Archive for 

Rational Mechanics and Analysis 130 (1995), 183-204. 
[Osh73] Stanley Osher, Initial-boundary value problems for hyperbolic systems in regions with corners. I, 

Trans. Amer. Math. Soc. 176 (1973), 141-165. 
[Osh74] , Initial-boundary value problems for hyperbolic systems in regions with corners. II, Trans. 

Amer. Math. Soc. 198 (1974), 155-175. 
[PT10] M. Petcu and R. Temam, The one dimensional shallow water equations with Dirichlet boundary 

conditions on the velocity, Discrete and Continuous Dynamical Systems Series S 4 (2010), no. 1, 

209-222. 

[PT11] , The one- dimensional shallow water equations with transparent boundary conditions, 

Math. Meth. Appl. Sci. (2011), DOI: 10.1002/mma.l482. 

[RTT08a] A. Rousseau, R. Temam, and J. Tribbia, Boundary value problems for the inviscid primitive 
equations in limited domains, in Computational Methods for the Atmosphere and the Oceans, 
Special Volume of the Handbook of Numerical Analysis, vol. XIV, Roger M. Temam and Joseph 
J. Tribbia, Guest Editors, P.G. Giarlet Editor, Elsevier, Amsterdam, 2008. 

[RTT08b] , The 3D Primitive Equations in the absence of viscosity: Boundary conditions and well- 

posedness in the linearized case, J. Math. Pures Appl. 89 (2008), 297-319. 

[Rud91] W. Rudin, Functional Analysis, Second ed., International Series in Pure and Applied Mathemat- 
ics, McGraw-Hill Inc., New York, 1991. 

[SLTT] M.-C. Shiue, J. Laminie, R. Temam, and J. Tribbia, Boundary value problems for the shallow 
water equations with topography, Journal of Geophysical Research-Oceans 116 (2011), C02015. 

[Tan78] Masaru Taniguchi, Mixed problem for wave equation in the domain with a corner, Funkcialaj 
Ekvacioj 21 (1978), 249-259. 

[TT03] R. Temam and J. Tribbia, Open boundary conditions for the primitive and Boussinesq equations, 
J. Atmospheric Sci. 60 (21) (2003), 2647-2660. 

[WPT97] T. Warner, R. Peterson, and R. Treadon, A tutorial on lateral boundary conditions as a basic and 
potentially serious limitation to regional numerical weather prediction, Bull. Amer. Meteor. Soc. 
(1997), 2599-2617. 

[Yos80] K. Yosida, Functional Analysis, Sixth ed., Springer- Verlag, Berlin, 1980. 

The Institute for Scientific Computing and Applied Mathematics, Indiana University, 831 
East Third Street, Rawles Hall, Bloomington, Indiana 47405, U.S.A. 
E-mail address: AH: aimhuang@indiana. edu 
E-mail address: RT:temam@indiana. edu 



